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Pulsation at the tip of the first giant branch? 
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ABSTRACT 

The first results of our ongoing near-infrared (NIR) survey of the variable red giants 
in the Large Magellanic Cloud, using the Infrared Survey Facility (IRSF) and the 
SIRIUS infrared camera, are presented. Many very red stars were detected and found 
that most of them are variables. In the observed colour-magnitude diagram (J — K, K) 
and the stellar K magnitude distribution, the tip of the first giant branch (TRGB), 
where helium burning in the core starts, is clearly seen. Apart from the genuine AGB 
variables, we found many variable stars at luminosities around the TRGB. From this 
result, we infer that a substantial fraction of them are RGB variables. 
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1 INTRODUCTION 

After leaving the main sequence, stars of low to intermediate 
mass (0.8M© < M < 8Mq) become red giants and climb the 
first giant branch (RGB). At the tip of the RGB (TRGB), 
they experience the so-called helium flash, followed by a sud- 
den decrease in luminosity. Then, the stars become red gi- 
ants second time and ascend the asymptotic giant branch 
(AGB), increasing the scale of mass- loss and pulsation. On 
the upper AGB they become optically obscured because of 
the circumstellar dust shell formed by heav y mass- loss, mak- 
ing t hem accessible only in the infrared (Iben & Renzini 



1983) 



Recently, large and homogeneous NIR photometry data 
for the Large Magellanic Cloud (LMC) has becom e available 
owing to la rge-scale NIR surveys such as 2MA SS ( |5krutski~ 
et al] [19971 ) and DENIS ( |Epchtein emT]|l997| ). However, as 
they are single-epoch surveys, they do not provide variabil- 
ity information. Variability information is necessary to get 
accurate mean magnitudes that are also indispensable for a 
proper understanding of the red giant evolution. 

Many optically- visible variables were found in the LMC, 
as the by-products of the gravitational lensing experiments 
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like MA CHO (IWoodl BOOOl) MOA flNoda et alj |2002| ) and 
OGLE ( Udalski et al.]p^999|) . However, most, if not all, of 
high mass-losing variables might have been missed by these 
projects, because they are not readily detectable in the opti- 
cal wave bands due to their surrounding circumstellar dust 
shells. Typical examples of such dust-enshro uded variable s 
are OH/IR stars and dusty carbon stars (e.g., Habing 1996| ). 
As these dust-enshrouded variables are thought to be at the 
latest stage of AGB evolution (e.g., Habing 1996), the under- 
standing of final evolution of low to intermediate mass star 
is greatly helped if their accurate quantities such as their 
number density and other physical parameters (mean mag- 
nitudes and periods etc.) would be obtained through deeper 
and larger survey. 

Many low- amplitude variables were newly discovered by 
the aforementioned optical surveys. Th e vast major i ty are 



at the luminosities below the TRGB. Alves et al 



(1998 



suggested that they are AGB stars in the early evolution- 
ary phase prior to entering the thermal-pulsing phase of the 
AGB. However, it is not clear yet if they really are or not. 

In this letter, we present some first results of our NIR 
monitoring survey which was started in December 2000, 
using the InfraRed Survey Facility (IRSF) at Sutherland, 
South African Astronomical Observatory. The limiting mag- 
nitudes (10cr) of this survey are 17.4, 17.0 and 15.8 at J, H 
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Table 1. Typical photometric errors 



magnitude 


J 


H 


K 


10.0 


0.018 


0.018 


0.018 


11.0 


0.018 


0.018 


0.018 


12.0 


0.018 


0.018 


0.018 


13.0 


0.018 


0.018 


0.020 


14.0 


0.022 


0.021 


0.031 


15.0 


0.030 


0.029 


0.049 


16.0 


0.048 


0.048 






Figure 1. A schematic illustration showing the survey region. 
The background image is a visible one and the coordinates are 
based on ICRS 2000 system. Each mesh has a side of 20 arcmin 
and consists of 9 sub-meshes, which correspond to the field of 
view of the SIRIUS. 



and K, respectively. At these limits, we can cover all of the 
AGB stars except extremely red ones (i.e., J — K > 5.0, such 



as th ose found in NGC419 and NGC1978 by Tanabe et al 



1998) and partway down the RGB populations in the LMC. 



We have been monitoring the total area of 3 square degrees 
along the LMC bar (see Figure. |l|), sufficiently large to do 
statistical analysis and to make complete catalog of variable 
red giants. 



2 OBSERVATIONS AND REDUCTIONS 

The IRSF consists of a specially constructed 1.4m alt- 
azimuth telescope to which is attached a large-format 
3-channel infrared camera, SIRIUS (Simultaneous three- 
colour InfraRed Imager for Unbiased Surveys). The SIRIUS 
can observe the sky in the three wave bands J(1.25/im), 
if (1.63/im) and i^ s (2.14/im) simultaneously and has a field 
of view of 7.8 arcmin square with a scale o f 0.45 arcsec/pixel. 
Detail s of the instrument are fou nd in Nagashima et al 
(|l999|) and |Nagayama et al.| (|2002|). 



2.1 Photometry 

Standard image reduction (i.e., dark subtraction, flat-field 
correction, sky subtraction) was made by using the SIR- 
IUS pipeline (Nakajima, private communication). All images 
were reduced in the same way and a single image comprises 
10 dithered 5 sec exposures. Photometry was performed on 
the dit hered images with DoPHOT in the fixed position 
mode ( Schechter, Mateo fc Saha| 1993 ). 

At first, the instrumental magnitudes were transformed 
to the LCO system ones bas ed on observations o f a few dozen 
stars from Persson et al. (Persson et al. 1998| ) . Then, the 



LCO system magnitudes were further transformed to the 
CIT/CTIO system ones through the equations in Persson et 
al. We estimate the uncertainties in the series of these system 
transformations should not exceed 0.05 mag. The photom- 
etry has a signal-to- noise ratio (S/N) of about 20 at J, H 
and K magnitudes of at 16.1, 16.1 and 15.0, respectively. 
Typical photometric errors are summarized in Table [lj 

In the present study, the data are not dereddened. In- 
stead, we show the direction and magnitude of the reddening 
vector in each diagram, based on the relations in Koornneel 



( |l982| ). 



2.2 Detections of variable stars 

We used Image Subtraction Method ( ISM) package 



ISIS.V2.1 (|Alard fe Lupton|[l998| ; |Alard||200C| ) to detect vari- 
able stars. This method can find variable stars even in very 
crowded fields such as LMC bar centre, and its high effi- 
ciency in detecting variables is well k nown from the pre- 
vious/ongoi ng: surveys (MA CHO, e.g., Alcock et al.| 2001 ; 
MOA, e.g., [Bond et alj|200l[ OGLE, e.g., |Alard||l999[ ; 
In this study, 



variable stars with light 
mag or more would all have been detected 
ten thousand variables have been detected in our surveying 
area, and their JHK light curves have been accumulated. 



etc.). 

variations of 0.03 
So far, nearly 



3 RESULTS 

In the survey area, 717,239 sources were detected by 
DoPHOT. Because of the scan overlaps, there should be 
many multiple entries of a single source in this original sam- 
ple. We eliminated such multiple entries based on (1) spatial 
proximity (|Ar| < 1.5 ) and (2) difference in photometry 
(\Am j,HandK\ ^ 0.1 mag). It leaves 562,429 sources. Fur- 
ther, requiring detections at all wave bands with S/N^20 
in K leaves 226,793 sources. Among the 226,793 sources, 
DoPHOT classified 184,677 sources as 'perfect' stars in all 
three wave bands (i.e., object type = 1 in J, H and K). 
Hereafter, we refer these 184,677 'perfect' stars as the main 
sample. 

Among the main sample, 5,188 stars were found to be 
variables. The full discussions of their pulsation properties, 
such as pulsation period, amplitude, variability type and 
mean magnitude etc. will be made in the future paper. In 
this letter, we will concentrate only on the information that 
a star is variable or not. 
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Figure 2. Left panel: Colour-magnitude and colour-colour diagram of the 184,677 main sample stars. Right panel: The same as the 
left panel, but of the 5,188 variables found in the main sample. Each plane is binned by 0.025 mag and the number density levels are 
logarithmic (see the wedge). The Reddening vector is drawn for Eb-v — 1-0. A dashed diagonal line indicates the lower detection limit. 



4 DISCUSSION 

In Figure |2| left panel, we show the colour-magnitude dia- 
gram (J — K, K) and colour-colour diagram (H — K, J — H) 
of the main sample. In each diagram, numbers of stars in 
bins of area 0.025x0.025 mag 2 are computed and the fidu- 
cial colour is applied according to the number density of 
stars in each position (see the annotated colour wedge). Be- 
cause the plotted magnitudes are single-epoch ones, these 
are affected by the scatter due to the light variation of each 
source. As one can see, the lower right area and the upper 
part (K < 9.5) of the colour-magnitude diagram are sparsely 



populated. These are attributed to the J limiting magnitude 
and K band saturation of the SIRIUS, respectively. 



Nikolaev & Weinberg ( pOOq ) analyzed the 2MASS LMC 
colour- magnitude diagram and identified the major stellar 
populations. Compared to their colour-magnitude diagram, 
ours has a narrow and sharp shape, suggesting a smaller 
Galactic foreground content in our sample. As this indicates, 
we can reasonably ignore the contribution of Galactic fore- 
ground stars here, because our survey region is centred on 
the LMC bar (see Figure. |l|), where foreground contribution 
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is negligibly small compared t o the number of the intrinsic 
LMC stars ( |Cioni et al.]|2000a| )- 

From the result of our variable star survey, now we know 
which star is a variable. Since neither 2MASS nor DENIS 
provide variability information, it will be interesting if we 
add it to the colour-magnitude and the colour-colour dia- 
gram. Figure |2| right panel is the same diagram as the left 
one, but it only contains the 5,188 variables extracted from 
the main sample. By comparing left and right panels, we can 
tell which part of the diagram the variable stars reside in. 



4.1 Obscured variables 

Stars in the dust-enshrouded phase are particularly inter- 
esting because of their importance in many areas, such as 
in the supply source of interstellar matter and as a key to 
answ er long sta nding problems, like the carbon star mys- 
tery (Iben 1981). Because of their redness, most, if not all, 



of these red variables have been possibly missed by the pre- 
vious optical surveys, such as MACHO, MO A and OGLE. 
Previous o bservations o f dust-enshrouded variables in the 
LMC (e.g., IZiilstra et al 



1996; Reid 



1991 



Wood et al 



1992 



van Loon et al.| |1998| ) were biased toward bright ones such 
as IRAS sources, and were carried out only for a few tens of 
stars. Therefore, it will be invaluable if their pulsation pe- 
riod, amplitude, and mean magnitudes would be known by 
our survey. We have to defer the full discussion of the prop- 
erties of the obscured variables pending further variability 
observations, but it is of interest to show some preliminary 
results. 

The Midcourse Space Experiment (MSX) satellite sur- 
veyed the whole LMC in the mid -infrared and many dust- 
enshrouded stars were dete cted (Egan, van Dyk & Price 
200l|; |Wood fe Cohen] |200l[). These red stars are thought 



to be obscured by the dusty circ umstellar envelope , which 
was developed by their mass-loss (Tanabe et al. 1997). In the 
main sample, there are 254 very red (J — K > 2.4) stars. 
Most of them might be i nclud ed in the MSX sources anal- 
ysed by |Wood fe Cohe^| fl200l| ). By comparing the left and 
right panel of Figure. ^| one can easily notice the fact that 
most of such mass- losing stars are variables (194 stars out 
of 254 stars, 76.4%). This observational result corroborates 
the idea that the ma ss- loss and variab i lity a re connected. 

Some stars from Hughes & Wood ( 1990| ) have also very 
red colour and most of them are Carbon-rich. However, their 
sample only extends up to J — K < 3.0, and there is still 
a possibilit y that dusty Oxyg e n-rich variables l ike O H/IR 



stars (e.g., Wood et al. 1992 ; van Loon et al. 1998| ) may 
have redder colours. Many obscured variables (up to J 
K ro 5.0) were found by our survey, but for most of them, 
the spectral type remains unknown. Infrared spectroscopic 
follow-up work would be interesting. 



4.2 RGB variables? 

In spite of the difficulty of distinguishing between the AGB 
and the RGB stars on the colour-magnitude diagram, vari- 
abilit y is usually regarde d as a sign of AGB stars. How- 



ever, |Origlia et al. ( J2QQ2 ) recently found a significant scale 
of mass-loss among the first giant branch stars near the RGB 
tip. Since a strong mass- loss is believed to be associated with 
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Figure 3. Stellar K magnitude distributions differentiated by 0.1 
mag bins of main sample (upper panel) and of 5,769 variable stars 
extracted from it (lower panel). Note that the tip of the RGB is 
clearly seen. 



pulsation, it is natural to search for the light variation in the 
RGB stars. 

In the upper panel of Figure [| we show the differential 
stellar K magnitude distribution N(m) (dK =0.1 mag) of 
the main sample. The discontinuity in the N(m) is clearly 
seen around K « 12.2, which corresponds to the TRGB. 
Most of the stars lying above the TRGB are intermediate 
age AGB stars, while both RGB and AGB stars are present 
just under the TRGB with the ascendancy of the RGBs. 



Renzini 



(1992) suggested that in a composite population 
the N(m) must exhibit a drop by more than a factor of 4 at 
the luminosity corresponding to the TRGB, and the ratio 
of AGB to RGB stars near the tip of the RGB is about 
1/3 for the intermediate age stars. Our result confirms the 
predominancy of the intermediate population in the LMC. 

We show the K magnitude distribution for the 5,769 
variables in the lower panel of Figure ^. As is clearly seen 
in it, there are two peaks in the luminosity of variables; one 
consists of genuine AGB variables and the other consists 
of stars just below the TRGB discontinuity. The position 
of the dip between the two peaks coincides with the above 
mentioned discontinuity for the main sample. 



Alves et al. (1998) and Wood et al. (1999) discussed the 



MACHO red variables that lie below the tip of the first giant 
branch. They suggested that all red giant vari a bles a re most 
likely AGB stars. According to |Origlia et al. ( |2002| ), about 
15% of RGB stars near the TRGB have dust shells. Assum- 
ing all AGB stars vary and all of these 15% of RGB stars 
near the TRGB also vary, and using Renzini's AGB/RGB 
ratio of 1/3, the ratio of variable RGB stars to AGB stars 
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is at least 0.45. Although the second peak might be made of 
AGB variables with much fainter luminosities than the typ- 
ical AGB variables, the coincidence of the dip luminosity for 
the variables with the discontinuity luminosity for the main 
sample can be more naturally explained by assuming that 
the second peak is due to a contribution from RGB variables 
since there is no reason to assume, according to models, that 
early AGB variables accumulate at the TRGB ( |Cioni et al. 
2QQ0b[). 



5 SUMMARY 

In this letter, we presented the first results of our ongoing 
NIR survey of variable red giants in the LMC. Many ob- 
scured stars (J — K > 2.4) were detected and we found that 
most of them are variables, indicating that the mass-loss 
and variability are connected. The nature of the obscured 
variables are still unknown and further observational inves- 
tigations are necessary and in progress. We also found that 
luminosity distribution of variable stars has two peaks. One 
is well above the TRGB and the other is just around it. 
Because there is no obvious reason to assume that Early 
AGB variables "pile- up" around the TRGB, we attributed 
the second peak to a contribution from RGB variables. From 
this point of view, we propose the evolutionary scenario of 
red giants as: (1) pulsate at the tip of the first giant branch; 
(2) pulsation fades away with the helium flash; (3) pulsation 
recommences as they ascend the asymptotic giant branch. 
However, the discrimination between RGB and AGB is still 
a challenging subject, and more direct evidence for their 
classification is needed. 
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